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Abstract 

We consider the motion of planar phase-transition fronts in first-order phase 
transitions of the Universe. We find the steady state wall velocity as a function 
of a friction coefficient and thermodynamical parameters, taking into account the 
different hydrodynamic modes of propagation. We obtain analytical approximations 
for the velocity by using the thin wall approximation and the bag equation of state. 
We compare our results to those of numerical calculations and discuss the range of 
validity of the approximations. We analyze the structure of the stationary solutions. 
Multiple solutions may exist for a given set of parameters, even after discarding non- 
physical ones. We discuss which of these will be realized in the phase transition as 
the stationary wall velocity. Finally, we discuss on the saturation of the friction at 
ultra-relativistic velocities and the existence of runaway solutions. 

1 Introduction 

In a first-order cosmological phase transition, bubbles nucleate and expand, converting 
the high-temperature phase into the low-temperature one (see, e.g., [HOE])- As bubbles 
expand, latent heat is released at their boundaries. This energy raises the temperature 
and causes bulk motions of the plasma. The perturbations caused in the cosmic fluid by 
the nucleation and expansion of bubbles generate a departure from thermal equilibrium. 
This may give rise to a number of cosmic relics, such as a baryon number asymmetry 
jl], baryon inhomogeneities [5], magnetic fields [6], topological defects [TJ, or gravitational 
waves 

In general, the system can be described by a relativistic fluid and a scalar field <p a ^ 
finite temperature T [101 [HI [12] . The latter may be a Higgs field and acts as an order 
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parameter. At high temperatures the free energy jF((j),T) has a minimum 0+(T) (in 
general, + = 0) and at low temperatures a different minimum 0_(T). In a first-order 
phase transition, the two minima coexist in certain range of temperatures, separated by 
a barrier. In the high-T phase, the free energy density is given by F + (T) = J r ((j) + (T),T), 
whereas in the low-T phase, it is given by J\_(T) = J-"(0_(T), T). The critical temperature 
T c is that for which J-" + (T c ) = J-"_(T C ). The phase transition occurs when the temperature 
of the Universe reaches T c . At T = T c , though, the nucleation rate vanishes, and bubbles 
effectively begin to nucleate at some temperature T n below T c [21 IT3] . 

The nucleated bubbles expand due to the pressure difference between the two phases. 
In general, the bubble walls reach a terminal velocity v w due to the friction with the 
surrounding plasma. Recently, the hydrodynamics of the moving walls has received much 
attention (see, e.g., [2J [TBI El EH]) due to the interest in performing thorough calcu- 
lations of the wall velocity and the energy injected into bulk motions of the fluid. These 
quantities are relevant for the generation of baryons and gravitational waves. In Ref. 
|15j . the ultra-relativistic velocity regime was considered, and it was shown that a state 
of continuous acceleration of the bubble wall is possible. Such "runaway" solutions may 
play an important role in the generation of gravitational waves. 

A realistic evaluation of the cosmological consequences of a phase transition requires 
considering the dynamics as completely as possible. Following the development of a phase 
transition involves the calculation of several temperature-dependent quantities such as, 
e.g., the pressure of the two phases and the bubble nucleation rate. During the phase 
transition the temperature varies in time and space due to the adiabatic cooling of the 
Universe and the release of latent heat. As a consequence, one has to deal with a set 
of integro-differential equations. Some of the involved variables are very sensitive to 
approximations (for instance, the nucleation rate). In order to avoid large errors, it is 
convenient to resort to nontrivial numerical calculations for these quantities. On the other 
hand, finding analytical approximations for other variables (e.g., the bubble wall velocity) 
provides a way of reducing the computation time. Widely used simplifications include the 
thin wall approximation and the bag equation of state PHI EE E21 EH] ■ Even with these 
approximations, it is not always possible to obtain analytical results. 

A considerable simplification is achieved by considering planar walls. Analytic results 
for the planar case were found recently in Ref. [H] for the wall velocity and in Ref. [17] 
for the energy injected into the fluid. It is important to note that considering spherical 
bubbles is not necessarily a better approximation than considering planar walls. Although 
the spherical symmetry is a good approximation for the initial stages of bubble growth, 
some cosmologically interesting outcomes of the phase transition are produced when bub- 
bles collide and lose the spherical symmetry. Moreover, losing the spherical symmetry 
is a requirement, e.g., for the generation of turbulence or of gravitational waves. As an 
explicit example, the "envelope approximation" for the generation of gravitational waves 
in bubble collisions neglects the overlap regions of colliding bubbles and follows only the 
evolution of the uncollided bubble walls. For such a calculation, the approximation of 
treating the walls as planar is, in principle, as good as considering spherical bubbles (but 
less complicated). In general, one does not expect important differences (see Ref. [17] for 
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a comparison of different wall geometries). 

In this paper we investigate the propagation of a planar phase transition front in 
the plasma. We aim at finding analytical approximations for the stationary velocity of 
the bubble wall as a function of the friction and the thermodynamical parameters. The 
present work is a continuation of the investigations of Ref. [14] . where we considered 
planar walls propagating as weak deflagrations or weak detonations. Here we include 
into consideration the case of supersonic Jouguet deflagrations [11 J and the possibility 
that the walls run away (see Ref. [16J for a recent study for spherical-symmetry walls). 
We also discuss here whether the different solutions are physical or not, and which of 
them will be realized as final stationary states during the phase transition. We discuss on 
the validity of the analytical approximations. The approximations are better for weaker 
solutions than for those close to the Jouguet points. For comparison, we consider some 
cases previously studied with numerical calculations [TUl [EH [T2] . 

The paper is organized as follows. In section [2] we consider the equations for the 
profiles of the fluid and the bubble wall, including a phenomenological friction term. We 
study the thin wall limit. In section [3] we use the bag equation of state to obtain a set 
of analytic equations for the wall velocity. In sections H] and [5] we present our results 
for the stationary motion and compare them with those of the numerical works of Refs. 
Pill El [12] ■ In section H] we discuss the range of validity of the analytical approximations 
and in section [5] we analyze the dependence of the wall velocity on the thermodynamic 
parameters and the friction. In section [6] we consider a different phenomenological friction 
term, which takes into account the fact that, in some models, the friction force approaches 
a constant in the ultra-relativistic regime. Finally, we conclude in section [71 

2 Hydrodynamics and microphysics 

We shall consider the motion of bubble walls at a given temperature T n < T c , i.e., we 
shall regard the nucleation temperature T n as a free parameter. Therefore, we shall not be 
concerned here with the calculation of the amount of supercooling. All the thermodynamic 
quantities (energy density, pressure, temperature, etc.) are derived from the free energy 
density. We have two phases, characterized by the minima <fr±. Thus, the equation 
of state (EOS) is different in each phase. For instance, the energy density is given by 
p± (T) = J-±(T) — TJ r ' ± (T), where a prime indicates derivative with respect to T. The 
pressure is given by p±(T) = — J-±(T). The enthalpy density is given by w± = p±+p±, and 
the entropy density by s± = w±/T. The speed of sound is given by c 2 ±{T) = dp±/dp± = 
p'±(T) I p'±(T) . The latent heat is defined as L = p + (T c ) — p_ (T c ). 

In the regions separating the two phases the scalar field varies from + to 0_. 
These interfaces are the walls of expanding bubbles, i.e., the phase-transition fronts. In 
general, temperature gradients arise during the phase transition, and the temperature 
varies beyond the bubble walls. Thus, the system is characterized by the scalar field 
0(x, t), the temperature T (x, t) of the plasma (which we treat as a perfect relativistic 
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fluid), and the fluid velocity v (x, t). These variables are governed by the equations 



d^-T—uW + g^Fj+d^d^ = 0, (1) 
d^<P + — + f l Tf(<P/T)u»d^ = 0, (2) 



where VtF" = (7,7V) is the four- velocity field of the fluid, with 7 = 1/a/1 — v 2 . Equation 
(JTJ) follows from energy-momentum conservation, d^T flv = 0, whereas Eq. ([2]) is the 
equation of motion for 0, where we have introduced a general phenomenological damping 
term proportional to <9 M in order to account for the friction force acting on the scalar 
field. In section E] we shall consider a modification of this term which does not grow with 
7 for large 7. The function / and the dimensionless friction parameter fj can be derived 
by considering the microphysics in specific models [HI EHl l2Tj - 

We are interested in the steady-state motion of bubble walls. Therefore, we assume 
stationary profiles moving with the wall at constant velocity. We shall consider planar- 
symmetry fronts moving in the x direction. In the rest frame of the front, all time 
derivatives vanish and Eqs. ([HE]) become 

— T|^7 2 t> = constant, (3) 
-Tf 7 V - T + I (S) 2 = constant, (4) 
S - % - vTf (HT) f x iv = 0. (5) 

For a given model the static, 1-dimensional Eqs. (EHS]) or ^he dynamic, 4-dimensional 
equations ([I}[2]), can be integrated numerically, e.g., using a lattice [TUl dH [12] ■ However, 
it is very useful to assume that the interface is infinitely thin, thus eliminating the scalar 
field profile from hydrodynamics considerations. Assuming a thin wall is in general a good 
approximation, as the wall width is much smaller than the width of the fluid profiles. In- 
deed, the latter is determined by the dynamics of the phase transition and is roughly given 
by the time scale t. In contrast, the width of the wall, determined by the characteristic 
length of variation of 0, is roughly given by the scale T~ x . In general, t is many orders 
of magnitude larger than T _1 . As we shall see, taking the thin wall limit in Eqs. (J3JH]) is 
trivial, whereas doing so in Eq. ([5]) requires additional approximations. 

2.1 Hydrodynamics 

Equations ([3]) and (j3J) relate the fluid variables on each side of the wall (in the rest frame 
of the wall). We shall use a + sign for variables just in front of the wall and a — sign for 
variables just behind the wall. Since vanishes outside the interface, we have 

w_w_7 2 = u> + f + 7 2 , (6) 
W-V 2 _'J 2 L+P- = w + v 2 + ^i\+p + . (7) 
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These equations give v + as a function of t>__. The solutions have two branches (see Fig. [I]), 
called detonations and deflagrations. For detonations the incoming flow is faster than the 
outgoing flow (\v + \ > \v-\). The value of \v + \ is supersonic in all the range < \v-\ < 1, 
and has a minimum at the Jouguet point \v-\ = c_. The minimum value of \v + \ is the 
Jouguet velocity Vj et , with Vj et > c + . For deflagrations we have \v + \ < \v-\, and \v + \ has 
a maximum value Vj ci < c + at the Jouguet point \v-\ = c_. The hydrodynamical process 
is called weak if the velocities v + and t>_ are either both supersonic or both subsonic. 
Otherwise, the hydrodynamical process is called strong. 




Figure 1: v + vs i>_ for the bag EOS (c+ = c_ = 1/V3), for a+ = 0.1. The upper branch 
corresponds to detonations and the lower to deflagrations. 

We may also have discontinuities in the fluid profiles out of the bubble wall, which 
are called shock fronts. In the reference frame of these surfaces Eqs. (JMZ]) apply. The 
functions w(T) and p(T) are the same on both sides of the surface (since we have the 
same phase), but the temperature is discontinuous. 

Away from the phase transition front, the field is constant and Eq. ([1]) becomes 

d t ( W1 2 -p)+d x ( W1 2 v) = 0, (8) 
d t {w 1 2 v)+d x [{w 1 2 v 2 +p)] = 0. (9) 

Since there is no characteristic distance scale in Eqs. ([HUH]), it is usual to assume the 
similarity condition [22], namely, that w,p and v depend only on £ = x/t. Using the 
relation p' = c 2 p', where a prime indicates derivative with respect to £, we obtain the 
equation for the velocity profile (see e.g. [17] ) 



(10) 
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Therefore we have either v' = 0, which gives constant solutions 



v i0 = constant, 



(11) 



or (£ — v)/(l — £v) = ±c, which gives two additional solutions. Of these, only the one 
corresponding to the + sign will be realized for the boundary conditions of our problem, 



which corresponds to a rarefaction profile. The speed of sound c in principle depends 
on T and, thus, may depend on £. From Eqs. (JSJED we also obtain the equation for the 
enthalpy profile, 



which can be readily integrated to obtain the profiles of the thermodynamical variables. 
In particular, for v = constant we have T = constant. 

The fluid velocity and temperature profiles are constructed with the solutions (limi3j) . 
using the matching conditions (JMZD and appropriate boundary conditions. The usual 
boundary conditions consist of a vanishing fluid velocity far behind the moving wall (at 
the center of the bubble) and far in front of the wall, where information on the bubble 
has not arrived yet. The temperature far in front of the wall can be determined from 
the dynamics of the phase transition. We shall assume it is given by the nucleation 
temperature T n of the bubble. Three kinds of solutions are compatible with all these 
requirements (see [T7j for details): 

A detonation. The wall is supersonic and the fluid in front of it is unperturbed 
(see Fig. |2J). Therefore, the fluid velocity in front of the wall vanishes, and we have 
v w = —v+, with |u + | > Vj ct . It turns out that the fluid profile is only compatible with a 
weak detonation. Thus, the outgoing flow is also supersonic, |i>_| > c_. Behind the wall, 
the fluid velocity is a constant between a certain point £o and £w = v w . At £ = £o the 
fluid velocity matches the rarefaction solution (ITS]) , which vanishes at £ = c_. For £ < c_ 
we have v — 0. As we shall see, the fluid profile behind the wall does not play a role in 
the determination of the detonation wall velocity. 

A "traditional" deflagration. For this solution the fluid behind the wall is at rest, 
so v w = — V-. The deflagration could in principle be weak, Jouguet, or strong (however, 
the latter seems to be unstable [TOl [TT| [T2]). The fluid velocity in front of the wall is a 
constant up to a shock front £ s h, which moves supersonically. Beyond the shock, the fluid 
is still unperturbed. 

A supersonic deflagration. The hydrodynamical solution is a Jouguet deflagration, 
i.e., < \v-\ and \v-\ = c_. Hence, we have \v+\ = Vj ei < c + . In this case, neither 
the fluid behind the wall nor in front of it is at rest, and the Jouguet condition t>_ = — c_ 
replaces the condition u_ = — v w of the traditional deflagration. Since the wall moves at 
the speed of sound with respect to the fluid behind it, and the fluid also moves (dragged 
by the wall), the wall velocity is always supersonic. The fluid profile behind the wall is 



«rar(0 



l-£c' 



(12) 




(13) 
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and £ w . In front of the wall, the fluid 



W — V J 



dot 



given by the rarefaction solution (fT2|) between c_ 
velocity is a constant up to the shock front. 

Fixing the thermodynamical parameters, the latter solution exists for c_ < v v 
As the wall velocity of the supersonic deflagration approaches the limit £ w = c_, the 
rarefaction wave disappears and this solution matches the traditional deflagration. As 
£ w increases, the shock front and the phase-transition front become closer. In the limit 



,dct 



the shock wave disappears and the solution matches the detonation. 
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Figure 2: Some fluid profiles for the bag EOS with a n = 0.1 and a c = 0.06 {T n /T c ~ 0.88). 
Solid lines correspond to a weak deflagration with v w = 0.3, dashed lines to a Jouguet 
deflagration with v w = 0.65, and dashed-dotted lines to a weak detonation with v w = 0.9. 
Left: the fluid velocity. Right: the fluid temperature. The dotted line indicates the 
critical temperature. 



2.2 Microphysics 

In order to obtain analytical results for the wall velocity, we shall implement the thin wall 
approximation in the friction equation fl5]). Following Ref. [H], we multiply Eq. (pjj) by 
dcj)/dx, integrate across the wall, and use (dJ r /d(f))(d(f)/dx) = dT jdx — (dF / dT)(dT / dx) 
HO]. We obtain 

P+-P--J (-^) dT - j vTf(<j>/T) (<P' f vidx = 0, (14) 

where <$' = d<p/dx. We will make approximations for the integrals in Eq. (JHJ), so that we 
can express the result in terms of the values of the variables outside the wall. 
First, we write the friction term as 

- J iFfWT) (<pf V1 dx = r/ Kl7+ ^ M7 ~ (15) 

(notice that the fluid velocity v in the rest frame of the wall is negative). The friction 
coefficient rj is usually approximated by its value in the limit T + = T_ = T n , i.e., when 
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the effects of hydrodynamics are neglected and only microphysics is considered. This is 
the limit in which the fluid is left unperturbed by the wall, and is valid for very weak 
solutions, for which v + = v_ = —v w . According to Eqs. (ll4H15p . the wall velocity is given 
in this limit by 

p.(T n ) - p + (T n ) 

u w /ui|micro j \- LVJ / 

V 



with 1] given by 



V = f}T I f^m^'fdx. (17) 



We shall consider rj as a constant in Eq. (1151) . which can be calculated for a specific model 
using the approximation (TTTj) . Similarly, for the first integral in Eq. (TT4")) we shall use a 
linear approximation for the entropy density inside the wall, which gives 

^w+ +s ->< r +- r ->. fl8) 



or 

With these approximations, Eq. ffT4"]) becomes [H] 

P+-P-- (T + _ T - ] + l V (|w+l7 + + |w - |7 - ) = °- (19) 

This expression is valid for any specific choice of the effective potential, since the informa- 
tion on the EOS is encoded in the variables p and s. Furthermore, considering a different 
phenomenological damping only amounts to modifying the function of v + and w_ in the 
last term and is straightforward, as we shall see in section [6] with an example. Therefore, 
Eq. ffT9]) is readily applicable to a large variety of models. 

The approximation (fTT)) neglects the variation of the fluid velocity inside the wall, 
since it is obtained from Eq. ffT5]) by replacing vy with its mean value. We expect this 
approximation to be better for weak solutions. Jouguet and strong processes cause larger 
perturbations of the fluid, and will give larger deviations. Indeed, as can be seen in Fig. 
[lj for weak solutions the difference between v + and t>_ is maximal at the Jouguet point. 
Similarly, Eq. (fl8|) approximates the entropy density inside the wall by the mean value 
(s + + s_)/2. We expect this approximation to be good for weaker solutions and to fail 
perhaps for Jouguet or stronger solutions, for which the difference between + and — 
variables can be large. As we shall see, since only weak or, at most, Jouguet processes are 
actually realized in the phase transition, the strongest deviations from the correct results 
will indeed occur for velocities which are close to the Jouguet deflagration or detonation. 

The friction force was calculated for some specific models for the case of the electroweak 
phase transition [21] . In the thin wall approximation, the wall profile can be estimated 
by neglecting the temperature variation and the last term in Eq. Thus, we have 

d(f)/dx = -^2A.F(0,T), with 

A^(0,T) = ^(0,T)-^(O,T), (20) 
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and Eq. (ITT]) gives 

77 ~ 77T / /(0/T)V2AJ-(0,T)#. (21) 

Roughly, we have i] m fjTa, where a — J <p' 2 dx is the surface tension of the bubble wall. 
This corresponds to setting /(0/T) = 1. The function /, though, can make a quantitative 
difference (see, e.g., [HJ ED])- However, since we are not going to consider any particular 
model but rather take fj as a free parameter, a different / in Eq. (121 p only amounts to a 
redefinition of fj. 



3 The bag EOS 

Equations flS]), (jTJ) and (lT9"j) can be solved for the wall velocity once the equation of state 
of the system is known. It is convenient to use as an approximation the bag EOS, 



jr + (T) = - a+ T 4 /3 + e, T_ (T) = -a_T 4 /3. 



(22) 



This simplification allows to find analytical expressions for the solutions. In this model 
the latent heat is given by L = As and the speed of sound is a constant c = 1/ v3 in both 
phases. It is customary to express the results as functions of the variable a = ej (a + T 4 ). 
As discussed in Refs. [14, 20J, for applications it is convenient to use L instead of e. 
Therefore, we define the parameters 



L 



OLr 



a h 



For the bag EOS Eqs. 



and (J7J give 



4a+Tl ' 



a, 



(23) 



6v- ~ 2 



6v- ~ 2 



o o 1 

+ a% + 3«+ - 3 



(24) 



which is plotted in Fig. [ljfor a + = 0.1. The plus sign corresponds to detonations and the 
minus sign to deflagrations. The friction equation ( TT9|) can also be expressed in terms of 
v + , f_, and aq_, 



4f + t>_a + 

1 — 3l7 + f _ 



2 / s_ 

3 V 



+ ^(| W+ | 7+ + |^|7-) = 0, (25) 



with 



and 



7\ 



1 — a_ 



1 + f + f _ 
1/3 — 



1/4 



(26) 



The latter expressions depend on the ratio a_/a + , which introduces a dependence on the 
parameter a c , 



1 — 3a r 



(27) 
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From Eqs. ( 1241125"]) we can find the velocities v + and w_ as functions of rj/L, a + and 
a c . However, the variable a + can be eliminated, since the temperature T + in front of 
the wall is a function of the nucleation temperature T n . The relation between a + and a n 
depends on the type of hydrodynamic solution. For detonations the temperature T + is 
just given by T + = T n , and hence a + = a n . For deflagrations, T + is related to T n through 
the matching conditions at the shock discontinuity, 

1 Vl 3T„ 4 + Tj 
viv 2 = -, — = 4 4 , 28 

where t>i is the velocity of the outgoing flow in the reference frame of the shock, and v-i 
that of the incoming flow. In the rest frame of the bubble center, the fluid velocity in front 
of the shock vanishes (see Fig. [2]). Hence, the velocity of the shock is given by t> s h = — 1>2- 
In the shock-wave region the fluid velocity is a constant and, thus, can be obtained either 
from the velocity v + or from the velocity v\. This gives the equation 

v w ~ \v+\ _ V3(a n - a + ) 

1-KK \/( 3 «™ + a+) (3a+ + a n ) ' 

Fixing the friction and the critical temperature, we can use the above equations to 
obtain function of a n as follows. For detonations we have v + = —v w . We can 

eliminate i>_ from Eq. (|24p . and then obtain v w from Eq. (|2"5|) as a function of a + = a n . 
For traditional deflagrations, we have u_ = — v w , so we can eliminate v + using Eq. fl24|) 
and obtain v w from Eq. (1251) as a function of a + . Then, Eq. (l29|) . together with Eq. (1241) . 
can be used to obtain a + as a function of a n and v w . For Jouguet deflagrations, we have 
v_ = —l/y/3 fixed, so Eq. (|24|) alone gives f + as a function of a + , i.e., f + = t>5 cf («+)• 
Therefore, Eq. (1251) gives already the value of a + (as a function of the parameters r\ and 
a c ). In this case, Eq. (I25|) can be used to obtain the wall velocity as a function of a + and 



4 General structure of the solutions 

Solving Eqs. ( 121H21?]) just amounts to finding the roots of algebraic equations, thus avoiding 
time-consuming numerical calculations. This is valuable when considering the develop- 
ment of a phase transition. However, in the way from Eqs. (J3J1S]) to Eqs. ( 1241129]) we 
have made several approximations. Therefore, it is important to determine the range of 
validity of these results. 

It is useful to consider the structure of the stationary states in the T + T_-plane. Follow- 
ing [11 j . we consider, on the one hand, the solutions which satisfy the energy-momentum 
conservation and the friction equation, i.e., Eqs. fl3HS]), but for which the boundary con- 
ditions have not been imposed. For the bag EOS and our analytical approximations, this 
amounts to considering Eq. ( 12 5p . with v + and t>_ expressed in terms of T + and T_, 

2 = (1/3 -a + -r/3)(l + q + + r/3) 

(1 + a+-r)(r + 1/3 - a + ) { ' 
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(1/3 - a + - r/3)(r - a+ + 1/3) 
(1 + ct+ - r)(l + r/3 + ct+) 



(31) 



with 



r 



a_Tt_ 



(32) 



This gives curves of rj = constant in the T + T_-plane (see Fig. [3]). On the other hand, we 
consider the solutions which satisfy the energy-momentum conservation equations, Eqs. 
(J3JHJ), and the boundary conditions, but for which the friction equation ([5]) has not been 
imposed. This gives curves of fixed T n in the T + T_-plane. For detonations, these are 
curves of constant T + = T n . For traditional deflagrations, the curves are given by Eq. 
(I2"9~j) with v w = |t>_|, and with v± given by Eqs. (|30ti32|) . Jouguet solutions are given by the 
condition v 2 _ = 1/3. The values of T + and T_ are constrained by the conditions < v± < 1 
and by the production of entropy at the phase transition front, s_|t>_|7_ > s + |t> + |7 + . 

The result is shown in Fig. [3j The dark grey region is forbidden by kinematics and 
by the condition of non-negative entropy production. There are two allowed regions. The 
upper one corresponds to deflagrations and the lower one to detonations. Weak deflagra- 
tions (white zone) are separated from strong deflagrations (light grey zone) by a dashed 
line which indicates Jouguet processes. The deflagration region is delimited by the line of 
zero entropy production (curved part of the boundary) and by the line of t>_ = (straight 
part of the boundary). The lower white region, corresponding to weak detonations, is 
delimited by the Jouguet condition (dashed line), the zero entropy-production condition 
(curved boundary), and the condition v+ — 1 (upper straight boundary). Solid lines 
represent the solutions of constant 77. Dotted lines represent the solutions for a given 
temperature T n for traditional deflagrations and for detonations. For detonations, the 
dotted curves are horizontal lines due to the boundary condition T + = T n . For the tra- 
ditional deflagrations, the curves are almost vertical lines, indicating that, although the 
temperature T + can be quite higher than T n due to reheating in front of the wall, the 
temperature T_ inside the bubble is very close to T n [TT] . If we plotted the dotted curves 
also for supersonic Jouguet deflagrations, they would lie on the dashed line separating 
weak and strong deflagrations. 

The possible stationary states correspond to the intersections of solid and dotted curves 
in Fig. [3J Thus, we see that for some values of the parameters (e.g., for T n /T c = 0.95, 
fj = 0.04) we will have multiple solutions for the wall velocity [TUl dH [12] ■ Notice also 
that solid lines in the deflagration region approach asymptotically the kinematic boundary 
corresponding to V- — 0. As a consequence, for large enough fj or for T n close enough to 
T c , there will always be weak deflagration solutions. Similarly, solid lines in the detonation 
region approach asymptotically the boundary of v + = 1. Therefore, there will always be 
detonations for small enough friction or strong enough supercooling. For intermediate 
values of fj and T n /T c , it may happen that neither weak deflagrations nor detonations 
exist, as e.g. for the case fj = 0.1, T n /T c = 0.89 (whereas for T n /T c = 0.86 we have a 
detonation and for T n /T c = 0.92 we have a weak deflagration). Generally, in such a case 
there will exist supersonic Jouguet deflagrations. 
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Figure 3: The solutions in the T + T_-plane for a + = (n 2 /90)g* with g* = 51.25, L = 
0.1T C 4 , and rj = rjTa, with a = 0.1T C 3 . The region with the dark shade is forbidden by 
kinematics or by the non-negativity of entropy production. The upper allowed region 
corresponds to deflagrations, and the lower one to detonations. The Jouguet processes 
are indicated by dashed lines, and the region with a lighter shade corresponds to strong 
deflagrations. Dotted lines correspond to fixed values of T n /T c = 0.86, 0.89, 0.92, 0.95 and 
0.98 as indicated, and solid lines to r) fixed (see explanation in the text). The red line 
corresponds to r) = 0. 



In Refs. [TQl [TT] [12], the differential equations for <p(x), T(x), and v (x) were solved 
numerically using a grid. Furthermore, a ^-dependent, quartic potential was considered 
as an approximation for the free energy. Comparing the results of those numerical com- 
putations with our results, we can test to what extent the use of the bag EOS, the thin 
wall limit and the approximations (I15H19P are valid. Thus, in Fig. [3] we have consid- 
ered "QCD-type" parameters used in Refs. [10J [HJ [12] and we have set f(cf>/T) = 1 
in Eq. ( 12T|) . which gives a friction of the form rj = r)To. The value of a + is given by 
a + = (7T 2 /90)(?*, where the number of effective degrees of freedom in the high-temperature 
phase is = 51.25. The latent heat is given by L = 0.1T C 4 , and the surface tension is 
cr = 0.1T C 3 . The set of values for the friction and for the nucleation temperature are the 
same as in Fig. 1 of Ref . [UJ . 
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The dotted lines, as well as the limits of the allowed region, are qualitatively identical 
and quantitatively quite close to those of Ref. [TT]. This indicates that the bag model 
is a good approximation for the EOS, at least for QCD-type parameters. In the white 
region (i.e., that of weak solutions) the solid curves are quantitatively similar to those of 
Ref. [llj. Qualitatively, the curves deviate from those of Ref. [TT] as they approach the 
dashed lines (i.e., as the solutions approach a Jouguet process). This is most apparent 
for fj — > 0. In the plot of Ref. [TT], for small friction the solid curves are parallel to the 
border of the allowed region (in the top-right corner of our Fig. E}. This is because the 
entropy production is related to gradients of through the friction parameter fj. Indeed, 
combining Eqs. (J3J1S]), one obtains 



As a consequence, the solid curve for fj = should coincide with the limit between the 
white and grey regions, corresponding to zero entropy production. Instead of that, in our 
case this curve (red line) enters the allowed region for detonations and enters the forbidden 
region for deflagrations. This is an indication of the break-down of the approximation for 
the integral of the entropy density across the wall, Eq. ( ITS]) , as weak solutions approach 
a Jouguet process (see the discussion in subsection 12. 2p . However, as we shall see in the 
next section, in most cases our analytical approximations give values of the wall velocity 
which are qualitatively and quantitatively good, even for solutions near the Jouguet point. 

In Ref. [TT] , the solid lines in the T + T_-plane do not penetrate the strong deflagration 
region but stop at the Jouguet line. This means that the numerical code did not find 
any strong deflagrations. As explained in Ref. [TT] , this is due to the fact that no matter 
how well one tries to guess the correct solution, the guess does not relax to a strong 
deflagration but rather changes considerably to form a different type of solution. This 
is an indication that strong deflagrations are unstable. In Fig. [3] it is seen that we find 
strong deflagration solutions. We believe this is not a shortcoming of our approximations, 
but rather due to the fact that analytical equations allow to find any solution, even the 
unstable ones. In this region, though, the departure of our solid lines from the actual 
curves is probably large. In any case, being unstable, strong deflagrations are of little 
interest. 



Let us consider now the wall velocity as a function of the parameters. As can be seen in 
Fig. HI weak traditional deflagrations (red curves) always exist for large friction (fj > 1) 
or little supercooling (a n ~ a c ). For lower values of the friction or larger amounts of 
supercooling, the traditional deflagrations surpass the speed of sound, becoming strong 




(33) 



which yields 




(34) 



5 The wall velocity 
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deflagrations. Shortly after crossing the sound barrier, though, the red curves end due to 
the condition of non-negative entropy production. At this point we always have supersonic 
Jouguet deflagrations (black lines), which match the traditional deflagrations at v w = c. 
For small enough friction or strong enough supercooling we always have detonations (blue 
lines) . 




Figure 4: The wall velocity for = 51.25, L = 0.1T C 4 , and i] = fjTa, with a = 0.1T C 3 . 
Detonations are plotted in blue, Jouguet deflagrations are plotted in black, and traditional 
deflagrations are in red. The dotted lines correspond to the sound and Jouguet velocities. 
Left: the wall velocity as a function of a n , for fj = 0.05. The range of a n correspond to 
values of the temperature between T n = T c and T n ~ 0.8T C . Thus, the a n -axis begins at 
the value a n = a c ~ 4.45 x 10~ 3 . Right: the wall velocity as a function of the friction, 
for T n = 0.95T C (a n ~ 5.46 x 10~ 3 ). 

Notice that, for some values of r] and a n , there are more than one stationary state. This 
can be seen already in Fig. [31 since the solid and dotted curves intersect at several points. 
In Fig. H] we show all the solutions for a given friction and supercooling. In some ranges 
of the parameters we have, for instance, a deflagration and a detonation, or a deflagration 
and two detonations. This is in agreement with the results of Refs. [101 [111 EG] • We may 
also have two traditional deflagration solutions for a given set of parameters, as shown in 
the left panel of Fig. [5J One of these solutions is a strong deflagration. 

One could expect that it will be always possible to choose the solution in such a way 
that, as the parameters are varied, one can go continuously from a subsonic, weak defla- 
gration to a supersonic, Jouguet deflagration and then to a weak detonation. As seen in 
Figs. H] and one can indeed change continuously from a weak deflagration to a super- 
sonic Jouguet deflagration. However, the Jouguet deflagration velocity does not match 
the detonation velocity. The continuity of the deflagration solutions is a consequence of 
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Figure 5: The same as in Fig. H] (right), but for T n = 0.891T C . In the right panel only one 
of the multiple solutions has been chosen. The dashed line corresponds to the velocity of 
the shock front. 

the continuity of the profiles. Consider the supersonic Jouguet deflagration. As the veloc- 
ity approaches the speed of sound c, the rarefaction wave vanishes continuously. Both the 
height and the width of the rarefaction vanish, matching continuously the profile of the 
traditional deflagration. In contrast, as the velocity approaches the Jouguet detonation 
velocity v j et , the shock wave becomes thinner and higher. Thus, the fluid velocity in front 
of the wall is maximal for a Jouguet deflagration at v w = Vj et (while the width of the 
shock wave vanishes). On the contrary, for a detonation the fluid velocity always vanishes 
in front of the wall. Therefore, the fluid velocity is discontinuous as the solution changes 
from a Jouguet deflagration to a detonation. This originates a jump in the parameters 
a n and fj as functions of v w , as can be seen in Figs. H]and|5j 

Leaving aside the jump of the parameters at v w = Vj et , one would expect that the su- 
personic Jouguet deflagration will always fill the velocity gap between weak deflagrations 
and weak detonations (as, e.g., in the left panels of Figs. SJandE]). A fluid profile does 
exist for any value of the wall velocity in the range c < v w < Vj ct . The Jouguet deflagra- 
tion reaches the value v w — c for the same values of parameters as the weak deflagration, 
whereas the value v w = Vj et is reached for a lower friction or a stronger supercooling than 
the detonation. This behavior can be explained by the fact that, for deflagrations, the 
compression wave which propagates in front of the wall and reheats the fluid, causes a 
friction effect [J3] which adds to the microphysics. Technically, the friction effect for the 
deflagration arises as a consequence of the relation ( 129|) between the nucleation tempera- 
ture T n and the reheated value T+. For the detonation, instead, we have T + = T n . This 
effect is easily estimated for weak solutions [13] • In the small supercooling limit (i.e., for 
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weak deflagrations) we have v w = Ap(T n )/r] eS , where Ap is the pressure difference p_ — p + 
and the effective friction is given by i] eS = (w-/w + ) [77 + L 2 / (\/3ii;_)l , which doesn't van- 
ish for 77 = 0. In the ultrarelativistic limit, in contrast, we have weak detonations with 
v w = 1 — 5, where 5 is proportional to rf and vanishes for vanishing friction. For deflagra- 
tions with high wall velocities (v w > c) the fluid velocity and temperature profiles have 
a sharp peak in front of the wall, so the friction effect can be considerably large. Notice 
that, in the right panel of Fig. [4J the value v j et is never reached by Jouguet deflagrations, 
even in the limit 77 — 0. Thus, we see that in some cases the effective friction can prevent 
the supersonic deflagration to reach the velocity of the Jouguet detonation. 

Even in the cases in which the Jouguet deflagrations fill the whole range between c and 
Vj et , only one of the multiple solutions will be realized. As a consequence, there will always 
be a gap in the velocity, as shown in Fig. (right panel). It is important to determine 
which of the possible solutions will be realized in the phase transition as a stationary state. 
In Refs. [TU| [TTJ [T2], the evolution from a given initial configuration of the scalar field and 
the fluid variables (i.e., from an initial "bubble wall") was studied by numerically solving 
the time dependent partial differential equations. Thus, their dynamical code selected 
one of the possible final states. It was observed that the strong deflagrations, as well as 
the branch of detonations which are closer to the Jouguet point, are not realized in the 
evolution of the bubble wall. Moreover, even if given as initial conditions, these solutions 
transform into one of the other solutions, suggesting an instability [UJ. Furthermore, 
as seen in Figs. |4] and [5j these solutions have a non-physical behavior as functions of 
the parameters [14]. Therefore, they must be discarded. In contrast, weak deflagrations, 
Jouguet deflagrations, and the branch of weaker detonations are in general stable. For 
high friction or low supercooling we can only have weak deflagrations. As the parameters 
are varied and the weak deflagration reaches the velocity v w = c, we must change to a 
Jouguet deflagration. On the other hand, for v w > c we will have to choose between 
a deflagration and a detonation. Both are stable and both are reachable from different 
initial conditions in the dynamical code of Ref. [UJ. However, in the normal evolution of a 
bubble wall, the solution which is realized is the detonation^- Accordingly, we must always 
choose the detonation if possible; otherwise, the Jouguet deflagration; otherwise, the weak 
deflagration. The result is shown in the right panel of Fig. |5l where we reproduced the 
plot of the left panel, keeping only the selected solutions. 

Qualitatively, our results agree with the numerical solutions, even for spherical bubbles 
(cf. Fig. 3 of Ref. [TJ] ). The parameters we used correspond to Fig. 13 of Ref. [TU] . 
There, the friction parameter is T — 1/fj. For a better comparison we plot the velocity 
as a function of T in Fig. [6j For the range of friction considered in that figure, we 
have a difference of at most a 5% for weak deflagrations and less for detonations. For 
the weakest solutions the difference vanishes, as expected. Moreover, for the strongest 
solutions, i.e., those around the Jouguet point, the agreement is quite good, as the jump 

1 Interestingly, the solutions which are not realized are those which are closer to the speed of sound and 
have a sharp peak in temperature. Thus, the dynamical evolution selects the weaker stationary solution, 
i.e., the weak detonation. For the fastest deflagrations that are realized, the bubble wall first goes into a 
detonation configuration before settling into a deflagration [T2] . 
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in v w is approximately at the same place, T m 10. 
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Figure 6: The wall (solid) and shock (dashed) velocities as functions of T = 1/fj. The 
values of the parameters are as in Fig. 

One expects that for any set of values of the parameters there will exist at least one 
solution. Consider a fixed value of T n . For large friction there should always be weak 
deflagrations and for small friction there should always be weak detonations. In the inter- 
mediate range there should exist Jouguet deflagrations. Although this is in general true, 
we find some exceptions in extreme cases. For instance, for a large amount of supercooling, 
T n /T c = 0.6 (which is out of the range of Fig. |3]), we find that, as we decrease the friction, 
the Jouguet deflagration velocity reaches the value of the Jouguet detonation before the 
detonation solution appears (see Fig. [JJ left panel). As a consequence, there is a range 
of values of rj for which there is no solution. For even stronger phase transitions, we find 
that there may be no solutions for large values of the friction. This behavior was observed 
also in Ref. [16J. This is shown in the right panel of Fig. [TJ where we considered values 
of the parameters similar to those used in Ref. |16j . For larger amounts of supercooling, 
there may be no subsonic deflagrations at all, whereas supersonic deflagrations cease to 
exist at some maximum friction. 

We have checked that, for physical models, the parameters (particularly the amount 
of supercooling) hardly fall into the case of Fig. [71 In particular, we considered the elec- 
troweak phase transition for several extensions of the Standard Model [9]. The situation 
of the right panel never arose, and that of the left panel arose only in a few limiting 
cases which are quite unlikely. For instance, extra scalars with very strong couplings with 
the Higgs field may yield an exceedingly strongly first-order phase transition. Only for 
some particular sets of parameters, and for the highest values of the couplings of the 
extra bosons, we found the situation of no stationary solution. Such couplings are ex- 
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Figure 7: The wall velocity as a function of fj for a c = 4.45 x 10~ 3 and T n /T c = 0.6 (left 
panel), and for a c = 0.05 (a_/a + = 0.85) and a n = 0.5 (right panel). The dotted lines 
indicate the speed of sound and the Jouguet detonation velocity. 

treme in the sense that a little increase causes the phase transition to remain stuck in 
the high-temperature phase and the universe to enter a period of inflation. In spite of 
this, the consequences of such models may be interesting and deserve further considera- 
tion. The absence of stationary solutions in some ranges of parameters may be related to 
the existence of runaway solutions [T5], which correspond to the wall propagating ultra- 
relativistically, with the gamma factor 7 growing linearly with time. Next we discuss this 
possibility. 

6 Friction saturation and runaway solutions 

The equation for the friction, Eq. fl2]), is obtained from the equation of motion for the 
field (ft, 

where V is the zero-temperature effective potential and is the distribution function of 
particle species i. The latter can be decomposed into the equilibrium distribution function 
/f q and a deviation 5 fa. The equilibrium part of fa in Eq. fl35|) gives dVr/dcf), where Vr is 
the thermal part of the finite-temperature effective potential (i.e., J 7 = V + Vt). Together 
with the term dV/d<j) this gives dJ-/d(f), which is the second term in Eq. (J2J), whereas the 
deviation gives the friction term. Usually, the deviations from equilibrium are assumed to 
be small, corresponding to a non-relativistically moving wall. This gives a friction force 
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proportional to the velocity of the wall with respect to the plasma. Such microphysical 
calculations can be used to determine the value of the coefficient fj and the function f(<p/T) 
in Eqs. ([2]) and (J5J). In this section we shall set for simplicity f(<f)/T) = 1. Therefore, this 
procedure gives a friction force per unit surface area 

= rjTav = rjv, (36) 

as explained at the end of section [2J Here, v is the (negative) fluid velocity in the wall 
frame. This corresponds to a friction term f}Tvd x <fi in the field equation (jSj). The simplest 
relativistic generalization is the u^d^ip term in Eq. (J2j). 

Recently, the opposite limit was considered in Ref. [15]. The friction acting on 
the electroweak bubble wall was derived for a wall which is already propagating ultra- 
relativistically, with gamma factor 7 ~ 10 9 . Such a fast moving wall validates a number 
of approximations. The reflection coefficients are exponentially suppressed. In the frame 
of the wall, incoming particles have received no signal that the wall is approaching and 
are in equilibrium. Furthermore, interactions or scatterings between plasma particles can 
be neglected and the occupancies evolve undisturbed. Therefore, only the equilibrium 
occupancies in the symmetric phase are used in the calculation. The occupancies are 
assumed to be constant along a classical trajectory through the bubble wall. 

As a consequence, the resulting force on the wall does not have a velocity-dependent 
term. The net force per unit area on the wall is found to be given by 

j = v(<p + ) - v{4>-) - - m*(<f> + )} I (2 ^, + /£(?)• (3T) 

This means that the total force per unit area is given by the "pressure difference" 

P--P+ = -£-+£+, (38) 

where F(4>, T) is the mean field approximation to the effective potential, which is obtained 
by keeping only the quadratic terms in a Taylor expansion of Vr about 

T(cj>,T) = V{4>) + V T (<P + ) + V[m^) - m?(0+)] CA 1 



^1 .vw lVY ^ J dm? 



(39) 

<{>+ 



To determine whether or not the wall can run away, T must be used instead of J 7 , and 
the total force ( 138]) must be positive, i.e., if replacing Vr(4>) with its second-order Taylor 
approximation removes the minimum <p_ or raises it above the minimum <p + , the bubble 
wall cannot run away [15J. 

In particular, in Ref. [15] it is shown that the bubble wall never runs away in a 
"fluctuation induced" first-order phase transition, i.e., a phase transition which is first- 
order due to the thermal part of the potential (e.g. the MSSM). As a simple example, 
consider the high-temperature expansion 

VtW = E _ mm + 0(m% m 
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It is well known that the cubic term in ( HO"]) may cause a first-order phase transition. 
This term is not present in the mean field potential. If the first-order character of the 
phase transition is due only to this term, then in the mean field potential the "broken 
symmetry" minimum 0_ ^ raises above the "symmetric" minimum <p + = 0. In such a 
model the wall will reach a terminal velocity v w < 1. An example of a model which does 
not rely on the terms Tmf (<f>) to yield a first-order phase transition is a potential with 
tree-level cubic terms. This is possible, e.g. in extensions of the Standard Model with 
singlet scalar fields, as considered in Ref. [To] . 

In this regime the friction force does not depend on the velocity of the fluid rela- 
tive to the wall, whereas for v — > it is proportional to v. It would be important for 
applications to know the behavior of the friction for intermediate velocities. Unfortu- 
nately, this constitutes a nontrivial problem. A simple interpolation between the two 
regimes was considered in Ref. [IB]. The approximation consisted in the replacement 
vPd„,<t> -> vPd^ /s/l + (A^) 2 in Eq. with A M = (0,0,0,1) in the wall frame. 
This is equivalent to the replacement wy d x (j) — > v d x <ft. This modification does not alter 
the discussion of section [2j Furthermore, it simplifies the analytic equations of section [3] 
through the replacement \v + \ 7 + + 7_ — > \v + \ + |t>_| in Eq. ff25l) . 

This phenomenological approach for the friction changes significantly the behavior in 
the detonation regime. As an example, we show in Fig. [S] how the results of Fig. are 
modified. Notice that detonations exist only in a small interval of the friction parameter. 




Figure 8: The wall velocity for the parameters of Fig. El with a friction force given by rjv 
(solid) and rj^yv (dashed). The dashed line is the curve of Fig. [5j 

For higher values of the friction we have deflagrations, whereas for smaller values there 
is no stationary solution and the bubble wall runs away. These results are in accordance 
with those of Ref. |16J . In particular, we obtain detonations only in a narrow region in 
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the ^a^-plane. 

On the other hand, the modification of the friction term does not affect significantly 
the deflagration solutions, as expected. In the previous section we have found regions of 
parameters corresponding to deflagrations, where there is no stationary solution (see Fig. 
[7]). In particular, the situation of the left panel of Fig. [7] may arise in some very strongly 
first-order phase transitions. With the modified friction, this behavior remains (see Fig. 
[H]). This is due to the fact that the deflagration solutions are not significantly altered by 
the modification of the friction. The runaway solutions appear instead for smaller values 
of the friction, corresponding to detonations. Nevertheless, the no-solution region has 
shrunk with this approximation, and it is possible that for a better approximation it will 
not exist at all. 
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Figure 9: The wall velocity as a function of fj for a c = 4.45 x 10~ 3 and T n /T c = 0.6, with 
a friction force given by r\v (solid) and lyyv (dashed). 

We wish to stress that this approximation, which consists essentially in assuming a 
friction force of the form ([36]) for any wall velocity, may be too simplistic. Although it 
reproduces the saturation of the friction force at large jv, the friction force saturates too 
soon (at values 7 ~ 1) to its ultra-relativistic value (which is in principle valid for 7 ~ 10 9 ). 
A different variation of the friction may allow the existence of stationary solutions (i.e., 
detonations) for a wider range of parameters. In particular, it is not clear that the friction 
should be a monotonically growing function of v. Indeed, the assumptions that lead to the 
runaway solution are based on the fact that the fluid is rather unaffected by the passage 
of the wall, due to the high speed of the latter. Conversely, the fact that the wall is 
unaffected by the fluid allows it to accelerate. This is somewhat similar to what happens 
at the macroscopic level with the stationary hydrodynamical modes, namely, the fastest 
the detonation, the weaker the disturbance it causes on the fluid [HI [T7J. Moreover, 
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solutions with intermediate velocities (around the Jouguet point) cause the maximum 
disturbance. Weaker deflagrations have small velocities and cause little perturbations. 
Weaker detonations cause less perturbations of the fluid and can thus move much faster. 
We may expect a similar behavior at the microscopic level, i.e., that intermediate velocities 
will cause larger departures from the equilibrium distributions and, thus, a higher friction. 

In particular, it may happen that the real friction force is well approximated by rjjv 
up to relatively high values of jv, such that 7 ^> 1 but not yet as large as to fulfill the 
hypothesis of Ref. [15J. If this is the case, then the wall will end up moving with a 
terminal velocity with a moderate value of 7 and never reach the ultra-relativistic regime, 
in spite of the existence of runaway solutions. In such a case, the model F&/A = rjjv 
considered in the previous section would give a better approximation than rjv, although 
the latter gives the correct ultra-relativistic behavior. Furthermore, it seems that this 
approximation does not take into account the fact that, if the first-order nature of the 
phase transition is fluctuation induced, then the bubble wall should not run away, even 
if the friction coefficient is very small. In such a case we expect the dashed line in Fig. |8] 
to give the correct behavior. 

7 Conclusions 

We have investigated the steady state motion of phase-transition fronts in a cosmological 
first-order phase transition. Our main goal was to find analytical approximations for 
the wall velocity, taking into account the different possibilities for the hydrodynamical 
modes and fluid profiles. Therefore, we have considered the case of planar walls, which 
allow to obtain analytical approximations. In Ref. [H] we considered the cases of weak 
deflagrations preceded by a shock front and weak detonations followed by a rarefaction 
wave. Here, we have studied also the case of Jouguet deflagrations which have both 
shock and rarefaction waves and move supersonically. We have considered two different 
phenomenological models for the friction. One of them grows linearly with the relativistic 
velocity jv [10], and the other saturates for large [16]. The latter reproduces the 
behavior of the friction force in the ultra-relativistic limit and leads to runaway solutions 

mi. 

Our main result is a set of algebraic equations which allow to obtain, from the thermo- 
dynamic parameters and the friction coefficient, the value of the wall velocity which will 
be realized as the final stationary state. These analytical results rely on several approx- 
imations, such as the use of the bag equation of state and the thin wall approximation. 
Implementing the latter in the equation for the friction requires some ansatz for the vari- 
ation of the entropy density inside the wall and also for that of the fluid velocity. The 
approximation for the entropy density seems to be the roughest one, as we do not obtain 
the curve of zero entropy production in the limit of vanishing friction. By comparing with 
numerical lattice calculations [TUl HU [12] , we have checked that the strongest departure 
from the exact solution occurs for the strongest physical solutions, i.e., those around the 
Jouguet points (either for detonations or deflagrations). 
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For a friction of the form jv, our results are in good quantitative agreement with the 
cases of planar walls considered in Refs. [TUl [EE] and in good qualitative agreement with 
those of spherical bubbles considered in Refs. [TTJ H2] . For a friction which saturates for 
large 7, the wall velocity shows essentially the same behavior as in Ref. [16], which consid- 
ered spherical bubbles. We remark that this latter approximation, although reproducing 
the correct behaviors for v — >■ and for v — >■ 1, may still be too simplistic for intermediate 
velocities. In particular, the friction saturates to a constant value for relatively small 
velocities, i.e., for values of the gamma factor which are much smaller than those which 
justify the approximations that lead to the runaway solution [15J. As a consequence, the 
region of parameter space in which the bubble wall runs away may be largely overesti- 
mated. The problem of determining the behavior of the friction at intermediate velocities 
is a difficult one and deserves further investigation, since it has important implications 
for cosmology. 
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